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Executive  Summary 


The  fraction  of  the  ocean  surface  covered  by  whltecaps 
has  long  been  thought  to  be  some  monotonieally  in¬ 
creasing  function  of  the  prevailing  wind  velocity  at 
least  for  large  fetches.  In  order  to  determine  the 
extent  to  which  other  factors  such  as  air  column  sta¬ 
bility  or  water  mass  type  can  influence  the  areal 
coverage  of  whitecaps,  photographic  data  was  collected 
over  the  Loop  Current  from  a  NAVOCEANO  P-3  aircraft. 
The  variation  of  whitecap  coverage  along  a  line  of 
closely  spaced  (25  km)  stations  was  determined  and 
compared  with  other  aircraft  and  data  buoy  informa¬ 
tion.  The  data  indicates  that,  on  the  day  of  the 
flight,  white-capping  within  the  boundaries  of  the 
Loop  Current  depended  little  on  the  local  wind,  but 
demonstrated  a  noticeable  dependence  on  air  column 
stability.  The  strength  and  nature  of  this  dependence 
varied  over  mesoscale  distances.  These  results  indi¬ 
cate  that  microwave  radiometric  measurements  can  be 
sensitive  to  variables  other  than  surface  wind  since 
microwave  brightness  is  quite  sensitive  to  sea  foam. 
The  value  of  areal  whitecap  coverage  as  a  measurable 
geophysical  variable  is  also  examined. 
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Background 

Wave  breaking  and  whitecapping  Is,  at 
the  same  time,  one  of  the  more  In¬ 
spiring  and  one  of  the  most  Important 
dynamical  processes  at  work  In  the  sur¬ 
face  layers  of  the  ocean.  To  the  wave 
modeler,  breaking  represents  a  sink  of 
wave  energy;  to  the  acoustic  scientist, 
a  source  of  noise  energy;  to  the  mixed 
layer  modeler,  a  source  of  turbulence  In 
the  surface  layer;  and  to  the  remote 
sensor,  a  bright  source  of  microwaves. 
Yet,  for  all  the  Importance  of  wave 
breaking,  It  Is  little  understood  and  is 
seldom  measured  In  the  deep  ocean. 

part  of  the  difficulty  of  analyzing 
whltecaps  stems  from  their  Intermittent 
character  In  both  space  and  time.  Since 
the  areal  coverage  of  whltecaps  typical¬ 
ly  represents  only  about  1%  of  the  sea 
surface,  an  Image  of  the  sea  surface 
must  encompass  on  the  order  of  500  m2 
in  order  to  accumulate  statistical 
samples  of  reasonably  small  (<  50%  of 
the  mean)  standard  deviation.  Eulerlan 
measurements  suffer  from  a  similar  sam¬ 
pling  problem.  bonguet-Hlggins  and  Smith 
(1981)  reported  that  measurements  made 
from  a  buoy  during  MARSEN  Indicated  that 
on  the  average  only  one  wave  broke  at 
the  buoy  during  every  100  significant 
wave  periods. 

Photographic  records  of  Instantaneous 
areal  whitecap  coverage  are  the  easiest 
whitecap  data  to  obtain,  and  several 
investigators  (Monahan,  1969;  Monahan, 
1971;  Toba  and  Chaen,  1973;  and  Ross  and 
Cardone,  1974)  have  attempted  to  relate 
whitecap  coverage  from  photographs  to 
measured  wind  stress  and  air  column 
stahllllty. 

Most  of  this  data  show ;  considerable 
scatter  when  whitecap  coverage  is 
plotted  against  wind  stress  alone. 
However,  Wu  (1980)  was  able  to  reduce 


the  scatter  of  Manahan's  (1971)  data 
by  organizing  the  data  according  to  the 
stability  of  the  marine  boundary  layer. 
Apparently  the  air  column  stability 
affects  the  formation  or  lifetime  of 
the  bubbles  that  constitute  the  foam; 
however,  no  present  theory  can  explain 
this  apparent  relationship.  Another 
source  of  scatter  is  the  lack  of  a 
well-defined  whitecap  boundary,  as 
discussed  in  Blanchard  (1971)  and  as 
demonstrated  in  Appendix  A  of  this 
report . 

In  parallel  to  the  field  measurements  of 
foam  coverage,  several  Investigators 
have  made  laboratory  studies  of  bubble 
formation  and  stability  as  functions  of 
the  physics  and  chemistry  of  the  water. 
Miyake  and  Abe  (1948)  measured  the  ef¬ 
fects  of  salinity  and  temperature  and 
found  that  for  oceanic  salinities,  bub¬ 
ble  stability  was  relatively  Insensitive 
to  the  small  changes  In  salinity. 

Changes  In  temperature,  on  the  other 
hand,  had  a  marked  effect  on  the  per¬ 
sistence  of  bubbles.  In  particular,  a 
5°  increase  In  temperature  at  20°C  pro¬ 
duced  a  5%  decrease  in  bubble  Lifetime. 
Garrett  (1967)  studied  the  effects  of 
surfactant  concentrations  on  the  stabil¬ 
ity  of  bubbles  and  found  that  small 
amounts  of  surfactants  Increase  bubble 
stability,  but  concentrations  in  excess 
of  one  raonomolecular  thickness  have  the 
effect  of  shortening  bubble  lifetime. 
There  have  also  been  in  situ  measure¬ 
ments  of  bubble  size  distributions  as  a 
function  of  depth.  Wu  (1981)  presents  a 
summary  of  much  of  this  work. 

Theory 

Monahan  (1971)  has  proposed  a  func¬ 
tional  form  for  the  areal  coverage  of 
whltecaps, 

W  =  A0  R^  (1) 
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where  W  is  the  fractional  coverage;  A0 
is  the  initial  area  generated  by  a 
breaking  wave;  R  is  the  number  of  break¬ 
ing  waves  occurring  per  unit  area,  per 
unit  time;  and  x  is  the  mean  lifetime  of 
a  whltecap.  x  has  been  measured  (Monahan 
and  Zietlow,  1969)  to  be  1.54  seconds 
for  fresh  water  and  3.85  seconds  for 
salt  water.  Hence,  the  instantaneous 
whltecap  coverage  is  a  measure  of  the 
four-second  time  history  of  breaking  as 
well  as  the  instantaneous  rate  of 
breaking. 

Wu  (1979)  reasoned  that  W  should  be 
proportional  to  the  energy  dissipation 
rate  of  breaking  waves  and  arrived  at  a 
dependence  of  W  on  the  wind  speed: 

W  ~  ul  ~  UlO3,75  (2) 

where  u*  is  the  friction  velocity 
and  is  the  10  meter  wind  speed.  If 

the  dependence  of  W  on  u*  in  (2)  is 
to  hold,  then  the  dependence  of  W  on 
other  variables  such  as  air  column 
stability  or  water  temperature  will 
likely  appear  as  an  additional  factor: 

W  ~  u^  f ( s ,  T,  etc..  .)  (3) 

where  s  is  some  measure  of  air  column 
stability  and  T  is  water  temperature. 
This  functional  form  is  suggested  by 
equation  (1)  if  one  assumes  that  the 
u^  dependence  resides  in  the  product 
AqR  and  the  other  independent  vari¬ 
ables  determine  x. 

In  this  paper  we  describe  an  experiment 
designed  to  examine  the  dependence  of 
whltecap  coverage  on  atmospheric 
boundary  layer  stability  over  the  ocean 
for  the  unstable  regime.  Whltecap  sta¬ 
tistics  are  measured  at  a  series  of  21 
consecutive,  collnear  stations,  each 
separated  by  23  kilometers  (km).  The 
line  of  stations  extends  across  the  Gulf 
of  Mexico  Loop  Current  along  a  line  of 
constant  latitude.  The  results  show 
that  for  slightly  unstable  boundary 
layers,  the  coverage  appears  to  decrease 
as  the  layer  becomes  more  unstable  in 
agreement  with  Monahan  (1971)  and  Wu 
(1980).  However,  for  extremely  unstable 
cases  the  whltecap  coverage  can  either 
Increase  or  decrease  with  stability, 


depending  on  the  locale.  The  whitecaps 
measured  were  generated  in  the  environ¬ 
ment  of  a  strong  current  and,  hence, 
may  not  be  representative  of  more 
benign  areas  of  the  ocean. 

The  Experiment 

An  instrumented  NAV0CEAN0  P-3  aircraft 
was  flown  over  the  Loop  Current  in  a 
westerly  direction  along  the  26°N 
latitudinal  line  as  shown  in  Figure  1. 

Da _a  were  collected  from  1 500S  to 
160n2  on  2/8/80.  The  following 
instruments  were  used  to  collect  data. 

PRT-5  Radiation  Thermometer  (SST) 
Thermal  Infrared  Scanner  (SST) 
Hasselblad  Camera  (Whitecaps) 

Laser  Surface  Profiler  (Sea  State) 
Wind  Anemometer  (Wind  speed/ 
direction) 

Air  Temperature  Sensor  (Air 
temperature) 

Expendable  BT’s  (Water  Temperature 
Profiles) 

The  flight  altitude  was  670  meters  (m) 
from  86°W  to  88°W  and  335  m  for  the 
remainder  of  the  track.  Aircraft  speed 
was  125  meters  per  second  (m/sec). 

Two  NOAA  Data  Buoy  Office  (NDBO)  data 
buoys  were  located  at  26°N,86°W  (2003) 
and  25. 9°N,89. 6°W  (2001),  respectively, 
as  shown  in  Figure  1.  The  buoys  supply 
hourly  reports  of  five  meter  height, 
three-rainute-averaged  wind  speed  and 
direction,  air  temperature,  water 
temperature,  and  barometric  pressure. 
These  two  buoys  were  the  principle 
source  of  surface  truth  data  for  this 
experiment . 

On  8  February  1980  the  Loop  Current 
flowed  northward  to  the  vicinity  of 
26 °N,  and  then  south  and  eastward  into 
the  Florida  Straits.  Figure  2  is  a  GOES 
Infrared  image  of  the  Loop  Current 
obtained  at  0600Z  on  7  February  1980 
that  shows  the  warm  core  of  the  Loop, 
the  cooler  slope  water  to  the  north  and 
the  east  and  a  warm  core  anticyclonic 
eddy  to  the  northwest.  The  sea  surface 
temperature  as  determined  by  the  air¬ 
borne  PRT-5,  and  verified  by  the  buoy 
reports,  is  shown  plotted  in  Figure  3. 
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Mexico  oircraft  track  with  surface  current  components 


surface  temperature,  whitecap  coverage,  surface  wind  stress,  and  RMS  waveheigh!  vs.  longitude 


The  sea  surface  temperature  difference 
between  the  Loop  and  the  Florida  shelf 
water  is  seen  to  be  approximately  5°C 
and  approximately  l.5°C  between  the  Loop 
and  the  eddy. 

The  Infrared  scanner  mounted  in  the 
aircraft  detected  two  fronts — the 
eastern  front  of  the  Loop  Current  and 
the  western  front  of  the  Loop  Current/ 
eastern  front  of  the  eddy  (Fig.  4). 


(1977).  A  reference  depth  of  350  m  was 
chosen,  and  the  resulting  surface  veloc¬ 
ity  components  are  shown  as  arrows 
along  the  aircraft  track  in  Figure  1. 
Similarly,  the  E-W  components  were 
derived  along  the  86°  meridian  and  are 
also  shown  in  Figure  1.  These  estimated 
values  for  the  surface  current  were 
later  used  to  adjust  the  air-sea  veloc¬ 
ity  shear  for  the  calculation  of 
surface  wind  stress. 


Figure  4  Infrared  scanner  images  of  loop  current  thermal  fronts 


These  frontal  boundaries  agree  with  ttie 
fronts  visible  on  the  satellite  image 
(Fig.  3). 

The  hydrographic  data  derived  from  AXBTs 
(Fig.  3)  shows  strong  frontal  character¬ 
istics  to  the  north  and  the  east  of  the 
Loop  Current  hut  weak,  gradients  to  the 
west.  The  N-S  components  of  surface 
velocity  along  the  track  were  estimated 
from  the  expendable  bathythermograph 
(XBT)  data  and  the  temperature-salinity 
(T-S)  characteristics  of  the  Florida 
Current  reported  by  Brooks  and  Siller 


Vhiteeaps  were  recorded  on  i  if  rare  i 
sensitive  film  (Kodak  2424)  and  were 
later  digitized  into  a  512x512  array  on 
an  image  processing  system.  Ten  images 
were  recorded  at  each  of  the  21  photo 
stations.  The  usable  field  of  view 
(FOV)  of  each  phot >3  was  370  m  x  560  m  up 
to  statLon  15  and  185  m  x  280  m  for  the 
remainder  of  the  track,  resulting  in  a 
pixel  size  oF  1  m  x  1  m.  The  image 
plane  of  the  camera  was  not  ..lotion  com¬ 
pensated  and,  lienee,  the  minimum  ground 
distance  resolvable  was  one  meter 


because  of  motion  blurring.  Twelve  of 
the  210  Images  are  shown  In  Figure  6. 
Streaks  running  north-south  were  evident 
In  a  few  of  the  photos.  One  such  photo 
Is  part  of  the  station  12  series  shown. 
Photos  containing  streaks  were  not  in¬ 
cluded  in  the  compilation  of  percentage 
whltecap  coverage.  The  changes  In 
whltecap  coverage  from  station  to  sta¬ 
tion  are  considerable  and,  as  can  he 
seen  from  Figure  6,  should  be  apparent 
to  the  most  casual  observer. 

The  digitized  Images  were  normalized  by 
accomplishing  an  intensity  stretching, 
area  by  area.  In  order  to  compensate  for 
uneven  lighting  distributions.  In  the 
event  tiiat  a  cloud  had  obviously  ob¬ 
scured  the  sun,  the  photo  was  discarded. 

Areas  of  each  image  that  appear  hright 
enough  to  be  whltecaps  are  "blotched"  or 
saved  on  a  one-dimensional  "graphics 
plane"  (1  x  512  array).  The  blotched 
areas  are  then  divided  by  the  total  area 
considered  to  obtain  the  fractional 
whltecap  coverage.  In  the  present  case 
all  pixels  with  an  intensity  value 
greater  than  160  on  a  0-256  scale  were 
blotched.  An  error  analysis  of  this 
somewhat  subjective  process  appears  in 
Appendix  A.  The  uncertainty  resulting 
from  such  an  analysis  Is  estimated  from 
Figure  A-3  to  be  approximately  +0.25% 
coverage.  Care  was  taken  to  assure  that 
digitizing  conditions  were  Identical  for 
all  210  Images. 

The  mean  percentage  whltecap  coverage  of 
each  group  of  Images  was  calculated,  and 
the  results  are  plotted  In  Figure  3  as  a 
function  of  longitude.  The  standard  de¬ 
viation  is  shown  by  the  length  of  the 
vertical  line  drawn  at  each  station. 

Many  images,  especially  near  the  western 
end,  were  affected  by  clouds  and  could 
not  be  used.  As  a  result  some  station 
means  are  based  on  as  few  as  three 
(e.g.,  station  19)  samples. 

Surface  Wind  Stress 

The  computer  program,  PR0PAR,  developed 
by  Cardone  and  reported  in  Cardone 
(1969)  was  used  to  Infer  the  wind  stress 
from  5  m  wind  velocities,  sea  surface 


temperature,  air  temperature,  barometric 
pressure,  and  an  estimated  relative 
humidity.  The  5  m  winds  present  at  sta¬ 
tions  between  the  data  buoys  were  ob¬ 
tained  by  fitting  the  aircraft-measured 
wind  speed  to  the  5  ra  winds  measured  at 
the  buoys.  This  amounts  to  assuming 
that  the  effect  of  the  higher  altitude 
is  to  scale  up  the  wind  velocities  by 
some  constant  factor  between  the  buoys. 
Sea  surface  temperature  was  measured 
from  the  aircraft  using  the  IR  radiom¬ 
eter.  Since  the  IR-derived  tempera¬ 
tures  at  the  buoys  were  approximately 
0.5°C  lower  than  the  In  situ  tempera¬ 
tures  measured  by  the  buoy  sensors,  they 
were  assumed  to  be  correct. 

Air  temperatures  measured  at  altitude 
were  reduced  to  the  surface  by  using  the 
buoy-measured  values  as  calibration 
points  In  the  same  manner  as  the  wind 
speed.  It  was  noted  that  the  tempera¬ 
ture  lapse  rate  was  higher  than  the 
adiabatic  lapse  rate  In  the  boundary 
layer  and  lower  at  altitude. 

Barometric  pressure  was  linearly  Inter¬ 
polated  between  the  buoys.  The  relative 
humidity  was  not  available  from  the 
buoys  at  the  time,  however,  buoys  pres¬ 
ently  instrumented  are  reporting 
between  50  and  100%  R.H.  A  value  of  70% 
was  assumed. 

Surface  wind  velocities  were  adjusted 
for  the  estimated  surface  current.  The 
surface  current  was  determined  through  a 
combination  of  the  surface  current 
components  derived  from  the  measured 
hydrography  and  the  pattern  of  surface 
Isotherms.  The  resulting  wind  stress  is 
shown  plotted  against  longitude  In 
Figure  3.  The  wind  stress  is  directed 
toward  the  west  and  northwest  and 
increases  raonotonically  from  east  to 
west,  primarily  as  a  result  of  the 
Increasing  wind  speed  (the  5  m  wind 
increases  from  8  ra/sec  at  the  east  end 
to  10  m/sec  at  the  west  end). 

Sea  State 

The  sea  state  was  characterized  by  an 
energetic  swell  propagating  from  the 
south  combined  with  shorter  wind  waves 
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Figure  6  Sample  images 


aligned  with  the  westerly  wind  direc¬ 
tion.  Evidence  of  both  wave  fields 
could  be  seen  in  a  Fourier-transformed 
photographic  image  of  the  sea  surface. 
Figure  7  shows  the  wave  height  power 
spectra  measured  at  the  two  NDBO  data 
buoys.  The  swell  period  of  6.25  seconds 
can  be  seen  in  both  spectra.  The  TV 
video  record  of  the  sea  surface  shows 
the  whltecaps  breaking  toward  the 
north,  indicating  that  the  northward- 
propagating  swell  was  the  dominant  con¬ 
tributor  to  the  observed  foam.  Since  the 
wind  turned  from  westerly  at  the  east 
end  of  the  flight  line  to  northwesterly 
at  the  west  end,  the  increase  in  swell 
amplitude  and  resulting  whitecapping  to¬ 
ward  the  west  can  be  explained  by  the 
interaction  of  the  northward-propagating 
swell  with  the  more  favorable  north¬ 
westerly  wind  in  evidence  at  huoy  2001. 


Figure  7.  Wave  height  spectra  at  NOAA  data  buoys 
2001  and  2003 

The  laser  profiler  on  the  P-3  records 
wave  height  as  a  function  of  along-track 
distance.  Due  to  periodic  noise  spikes 
in  the  data,  the  typical  continuous  wave 
height  record  length  was  of  order  3  km. 
The  profiles  were  digitized  and  root- 
raean-square  wave  height  computed.  The 
results  are  shown  plotted  against  longi¬ 
tude  in  Figure  3.  No  systematic  in¬ 
crease  in  significant  wave  height  ap¬ 
pears  from  east  to  west;  however,  the 


laser  records  did  not  extend  far  enough 
to  the  west  to  reflect  the  corresponding 
increase  in  the  whitecap  coverage  at  the 
west  end  of  the  flight  line. 

The  whitecap,  friction  velocity,  and 
air-sea  temperature  difference  data  are 
compiled  in  Table  I. 

TABLE  I 

Whitecap  No.  of 


Stat Ion 

Coverage 

Standard 

Ftiotos 

1*3. 

(*) 

Deviation 

AT 

Analyzed 

1 

0.57 

0.34 

-1.25 

6 

2 

0.55 

0. 29 

-2.90 

10 

3 

0.35 

0.13 

-5.15 

8 

4 

0.37 

0.08 

-4.65 

10 

5 

0.48 

0.11 

-3.85 

6 

6 

0.  83 

0.23 

-4.30 

7 

7 

0. 74 

0.13 

-3.65 

10 

8 

0. 82 

0.19 

-3.40 

9 

9 

0.66 

0.18 

-3.65 

8 

10 

0.54 

0.07 

-3.25 

10 

11 

0. 59 

0.09 

-3.05 

9 

12 

0.30 

0.07 

-3.10 

9 

13 

0.46 

0.19 

-3.40 

7 

14 

0.  54 

0.08 

-2.80 

9 

15 

0.47 

0.12 

-3.30 

9 

16 

0.41 

0.  11 

-3.15 

8 

17 

0.  73 

0.22 

-2.20 

9 

18 

0.84 

0.  11 

-1.55 

4 

19 

1.89 

0. 50 

-1.35 

3 

20 

1.06 

0.22 

-1.65 

7 

21 

Discussion 

1.46 

0.42 

1 

• 

yjt 

O 

7 

The  whitecap  coverage  shown  in  Figure  3 
appears  to  be  a  smoothly  varying  func¬ 
tion  of  space,  at  least  in  the  region  of 
the  Loop  Current.  Since  the  standard 
deviation  of  each  estimate  is  well  below 
the  mean,  the  fluctuations  are  believed 
to  be  geophysical  rather  than  a  result 
of  statistical  uncertainties.  The  most 
striking  property  of  the  whitecap  cover¬ 
age  is  its  weak  dependence  on  the  sur¬ 
face  wind  stress.  Whereas  the  wind 
stress  gradually  increases  from  east  to 
west,  the  whitecap  coverage  exhibits  a 
somewhat  wind-independent  behavior. 

Also,  coverage  appears  to  be  generally 
lower  within  the  Loop  Current  than  to 
the  west.  A  dependence  of  whitecap 
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coverage  on  Influences  other  than  the 
wind  is  suggested.  A  marked  decrease  In 
whitecap  coverage  occurs  Ln  the  vicinity 
of  both  Loop  Current  fronts.  A  similar 
decrease  in  coverage  was  observed  by 
Ross  and  Cardone  (1974)  near  the  north 
"wall"  of  the  Gulf  Stream. 

Whitecap  coverage  is  plotted  with  the 
air-sea  temperature  difference  as  an 
indicator  of  air  column  stability  in 
Figure  8.  A  definite  correlation  appears 
to  exist.  In  order  to  remove  the  wind 
stress  dependency  of  the  whitecap  cover¬ 
age,  W,  and  to  make  the  dependence  of  W 
on  stability  more  explicit,  we  scale  the 
whitecap  coverage  by  u*  as  suggested 
by  equation  (3),  and  plot  this  scaled  W 
against  AT.  The  resulting  plot  is  shown 
in  Figure  9. 


If  the  track  is  subsectioned  into  groups 
of  contiguous  stations  -  1-5  (A),  6-11 
(B),  12-16  (C),  and  17-21  (D),  as  shown 
in  Figure  9  -  a  localized  pattern  of  de¬ 
pendence  on  stability  emerges.  Sections 
A  and  D,  and  perhaps  C,  demonstrate  the 
dependence  of  whitecaps  on  stability 
concluded  by  Wu  (1981):  that  fractional 
coverage  decreases  with  decreasing  sta¬ 
bility.  However,  the  strength  of  that 
dependence  is  a  function  of  location. 
Section  C,  which  covers  much  of  the  warm 
core  of  the  Loop,  exhibits  an  increase 
in  whitecaps  with  decreasing  stability. 
While  the  entire  data  set  lacks  any  kind 
of  organization  with  respect  to  stabili¬ 
ty,  the  data  appear  to  be  well  organized 
locally  (within  100  km).  This  mesoscale 
variability  indicates  that  other  fac¬ 
tors,  perhaps  chemical  in  nature,  may 


90°  89°  88°  87°  86°  85° 


Figure  8.  Coverages  vs.  air-sea  temperoture  difference 
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affect  the  lifetime  of  the  whltecap  column  stability  were  shown  to  be  influ- 

bubbles  and  whltecap  coverage.  ential  in  determining  bubble  lifetimes 

and,  hence,  the  fraction  of  ocean  cov- 

The  fractional  whltecap  coverage,  W,  is  ered  by  foam  at  any  instant  of  time, 
shown  plotted  against  mixed  layer  depth  The  usefulness  of  such  information  de¬ 
in  Figure  10.  A  rough  correlation  exists,  pends  on  the  application.  For  example, 
showing  a  deepening  mixed  layer  with  aerosol  distribution  may  be  a  strong 

increased  whltecap  coverage  that  has  a  function  of  the  instantaneous  foam  cov- 

correlatlon  coefficient  of  about  0.5.  erage,  whereas  acoustic  noise  may  depend 

Again,  much  higher  correlations  can  be  more  on  the  number  of  whltecaps  gener- 

achteved  if  the  data  set  is  partitioned  ated  per  unit  time  per  unit  area. 

as  in  Figure  9.  Also,  the  locations  of 

the  two  minima  In  coverage  at  85°52'W  The  results  of  this  analysis  are  not  en- 

and  87a51'W  as  shown  in  Figure  3  corre-  couraging  for  the  prospects  of  abstract- 

late  approximately  with  the  locations  of  the  surface  wind  speed  from  the 

the  two  mixed  layer  depth  minima  shown  microwave  brightness  of  the  sea  surface 

In  Figure  5.  If,  as  assumed,  the  brightness  Is  deter¬ 

mined  In  large  part  (for  wind  speed  >10 

Conclusions  m/sec)  by  the  amount  of  foaming  taking 

place  on  the  sea  surface.  For  example, 

Spatial  fluctuations  in  whltecap  cover-  significantly  lower  microwave  brightness 
age  along  a  400  km  stretch  of  ocean  have  temperatures  might  be  observed  at  the 
been  shown  to  be  relatively  Insensitive  western  Loop  Current  boundary  than  would 
to  the  magnitude  of  the  local  surface  ke  observed  50  km  to  the  west.  This 

wind  stress.  Water  temperature  and  air  lower  temperature  would  be  Interpreted 
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as  Indicative  of  a  lower  surface  wind 
velocity  when,  In  fact,  the  reduced 
brightness  may  have  been  caused  by  a 
decrease  In  the  stability  of  the  bubbles 
caused  by,  say,  a  reduction  In  the  sur¬ 
factant  concentration  In  that  area.  The 
hoop  Current,  however,  may  represent  an 
area  of  anomalous  variability  In  the 
chemical  and  physical  properties  of  the 
sea  water.  These  conditions,  coupled 
with  the  hl-dlrect lonal  nature  of  the 
wave  field  and  the  limited  fetch  of  the 
wind,  combine  to  provide  what  may  likely 
he  an  atypical  situation. 

The  results  of  Appendix  A  have  shown 
that  whitecaps  are  not  the  well-defined, 
sharp-edged  features  they  appear  to  be. 
This  lack  of  whitecap  definition  limits 
the  utility  of  areal  whitecap  coverage 
as  an  independent  variable  In  functional 
expressions  for  wind  speed,  wind  stress, 
significant  wave  height,  etc.  Further 
research  may  show  that  some  Integral 
property  of  the  reflectance  of  foam  will 
prove  to  relate  more  closely  to  physical 
foam  parameters  such  as  bubble  spectra, 
bubble  concentration,  and  foam  depth 
and,  thus,  relate  more  closely  to  wave 
dissipation  rates,  aerosol  production, 
wind  speed,  and  microwave  brightness. 
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Appendix  A.  Errors  in  the  Estimation  of 
Areal  Whitecap  Coverage 
Using  Photographs 


A  cursory  glance  at  the  photos  displayed 
in  Figure  6  reveals  that  the  reflection 
of  light  from  foam  is,  hy  far,  the 
strongest  signal  present  in  a  photo¬ 
graphic  image  of  the  sea  surface  on  a 
sunny  to  partly  cloudy  day.  It  would 
seem  an  easy  task  (at  least  if  one  had 
access  to  an  imaging  computer)  to  meas¬ 
ure  the  area  of  the  photo  that  stands 
out  from  the  background  as  whltecaps  and 
divide  by  the  total  area  to  give  the 
fractional  whitecap  coverage.  The  in¬ 
tensity  of  the  image  of  each  whitecap  is 
not  raonotonic,  however,  but  varies  from 
a  maximum  value  near  the  core  of  the 
breaker,  tapering  on  every  side  down  to 
the  intensity  of  the  background.  In 
short,  the  boundary  of  the  image  of  the 
breaker  is  not  well-defined. 

Figure  A-l  is  a  magnified  portion  of  one 
of  the  images  used  in  the  foregoing 
whitecap  analysis.  The  whitecap  in  the 
left-center  of  the  picture  has  been 
digitized  and  the  intensities  contoured 
in  Figure  A-2.  While  the  center  of  the 
whitecap  has  limited  dynamic  range,  the 
edges  taper  gradually  in  Intensity  from 
the  bright  core  to  the  darker  back¬ 
ground.  The  cause  of  this  variation  is 
either  a  decrease  in  the  reflectance  of 
the  whitecap  as  the  boundaries  are 
approached  or  is  an  artifact  of  the 
exposure  process  whereby  the  photons 
incident  on  the  emulsion  are  of  such  an 
Intensity  that  scattering  processes 
within  the  emulsion  cause  a  lateral  dif¬ 
fusion  of  light.  The  degree  of  blooming 
Is  difficult  to  determine  a  posteriori, 
since  no  control  targets  were  Imaged 
during  the  experiment;  however,  Kodak 
specif ications  for  2424  film  indicate 
that,  for  contrasts  of  30  dB,  a  resolu¬ 
tion  of  80  tines  per  millimeter  (mra) 
should  be  obtainable.  As  seen  in  Figure 
A-2,  the  lateral  extent  of  the  whitecap 


image  on  the  camera  backplane  is  1.2  mm 
and,  hence,  descrlbable  by  some  100 
vertical  lines.  This  resolution  would 
certainly  be  sufficient  to  resolve  the 
observed  reflectance  variations  across 
the  whitecap  observed  In  Figure  A-2,  so 
we  must  conclude  that  the  measured  con¬ 
tours  are  representative  of  actual  re¬ 
flectance  variations  over  the  surface  of 
the  whitecap. 

Whatever  the  source  of  the  resulting 
intensity  variations,  the  problem  of  es¬ 
timating  the  boundaries  remains.  Several 
minimum  Intensity  criteria  for  defining 
a  whitecap  were  applied  to  a  sample 
image  (not  the  image  of  Figure  A-l),  and 
the  resulting  fractional  coverage  was 
plotted  against  the  minimum  brightness 
criteria  for  a  whitecap  boundary  (Fig. 
A-3).  The  minimum  brightness  is  selected 
from  a  full  dynamic  range  of  0-256.  The 
change  in  estimated  coverage  is  shown  to 
he  quite  large  for  a  corresponding 
change  in  threshold  brightness;  for 
example,  reducing  the  minimum  brightness 
criteria  from  160  to  155  produces  a 
two-fold  Increase  in  apparent  whitecap 
coverage . 

Since  the  threshold  value  Is  arrived  at 
quite  arbitrarily  (the  computer  analyst 
sets  the  level  at  whatever  value  seems 
to  best  define  the  whltecaps),  whitecap 
coverage  estimates  must  also  be  subject 
to  a  certain  amount  of  subjectiveness. 
Figures  A-4  and  A-5  demonstrate  the  ef¬ 
fect  on  the  blotched  areas  of  Increasing 
the  threshold  from  160  to  165.  These 
figures  are  enlarged  whitecap  Images 
overlaid  with  the  blotched  areas  for  the 
two  cases.  As  the  threshold  Is  raised 
In  Figure  A-5 to  165,  several  of  the 
pixels  from  Figure  A-4  no  longer  meet 
the  minimum  brightness  criteria  and  are 
rejected,  with  a  corresponding  drop  in 
the  blotched  area,  and  thus  In  the 
calculated  areal  whitecap  coverage. 
Hence,  variations  in  the  way  in  which 
areal  whitecap  coverage  is  determined  by 
different  Investigators  or  hy  the  same 
investigator  for  different  data  sets  can 
produce  such  large  variations  In  re¬ 
ported  whitecap  coverage  that  Inter¬ 
comparisons  may  become  meaningless. 
Blanchard  (1971)  arrived  at  this  same 
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Cipriano,  R.  J.  and  D.  C.  Blanchard 
(1981).  Bubble  and  Aerosol  Spectra  Pro¬ 
duced  by  a  Laboratory  Breaking  Wave.  J. 
Geophys.  Res., (in  press). 

Monahan,  R.  C.  and  C.  Zietlow  (1969). 
Laboratory  Comparisons  of  Fresh-Water  and 
Salt-water  Whitecaps.  J.  Geophys.  Res., 
v.  74,  p.  6961-6966. 


The  source  of  the  Intensity  variations 
within  each  whltecap  may  be  due  to  vari¬ 
ations  In  bubble  concentration  or  other 
properties  of  the  hubbies.  Cipriano  and 
Blanchard  (1981)  measured  bubble  spectra 
at  different  distances  from  the  base  of 
a  laboratory  waterfall  and  found  that 
the  concentration  of  bubbles  and  the 
fraction  of  bubbles  >300um  decreased 
with  distance  from  the  center  of  the 
foam.  The  observations  described  here 
indicate  that  the  reflectance  of  foam 
near  nadir  may  be  a  function  of  the 
shape  of  the  bubble  spectrum  and/or  the 
bubble  concentration. 

Obviously,  not  enough  is  known  about  the 
scattering  properties  of  sea  foam  and 
how  these  properties  are  related  to 
bubble  spectra,  bubble  concentration, 
foam  thickness,  foam  lifetime,  etc. 
Further  research  to  determine  the  rela¬ 
tionship  between  reflectance  and  the 
physical  properties  of  foam  will  deter¬ 
mine  which  quantities  other  than  areal 
coverage  -  Intensity  weighted  coverage, 
for  example,  might  be  more  closely  re¬ 
lated  to  energy  dissipation  or  microwave 
emissivity. 
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conclusion  after  comparing  his  results 
(Blanchard,  1963)  with  those  of  Monahan 
(1971). 

The  difficulty  In  judging  boundaries 
does  not,  however,  negate  the  value  of 
comparative  studies  on  the  same  data  set 
as  was  done  in  this  report.  That  is,  If 
the  criteria  for  establishing  foam 
boundaries  Is  kept  constant  together 
with  other  variables  such  as  digitizing 
illumination,  film  type,  filters, 
camera,  shutter  speed,  etc.,  meaningful 
comparisons  can  he  made. 


gure  A- 1  Magnified  sub-section  of  photo  image  of  vuhifecap  used  to  contour  reflected  light  intensity 


1  1 


UNCLASSIFIED 


fc  C  j  T  v  C  L  ASSi  ►  lC  ATiQN  Of  T  MlS  PAGE  fWhan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 


NORDA  Report  43 


4  TiTlE  'and  S.ibtltle\ 


?  GOVT  ACCESSION  NO. 


Measurements  of  Whitecap  Coverage  and  Surface 
Wubds  Over  the  Gulf  cf  Mexico  Loop  Current 


AUTHOR1'*; 


Peter  M.  Smith 


9  PERFORMING  ORGANIZATION  name  and  address 


Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  Mississippi  39529 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


HFC’PlENT'S  catalog  number 


Type  of  REPORT  A  PERIOD  COHERED 


PERFORMING  ORG  REPORT  NUMBER 


CONTRACT  OR  GRANT  NjMQER'%) 


10  PROGRAM  ELEMENT  PROJEC"  TASK 
AREA  &  WORK  UNIT  NjMBER$ 


t!  CONTROLLING  OFFICE  name  and  ADDRESS 


October  1981 


12  report  date 


Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  Mississippi  39529 


14  MONITORING  AGENCY  NAME  A  AOORESSf//  different  (torn  Controlling  Office)  15  SECURITY  CLASS  (of  thle  report) 


13  NUMBER  OF  P  AGES 

24 


Unclassif led 


15a  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  'of  this  Report) 


Distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  It  different  from  Report) 


19  KEY  WORDS  (Continue  on  reverse  aide  If  neceaamry  and  Idantlty  by  block  numbar) 


Remote  Sensing 
Whitecaps 

Air-Sea  Interaction 
Breaking  Waves 


Marine  Boundary  Layer 
Loop  Current 
Mixed  Layer 


20  ABSTRACT  (Continue  on  reverse  side  If  neceeeary  and  Idantlty  by  block  number) 

The  fraction  of  the  ocean  surface  covered  by  whitecaps  has  long  been  thought 
to  be  some  monotonically  increasing  function  of  the  prevailing  wind  velocity 
at  least  for  large  fetches.  In  order  to  determine  the  extent  to  which  other 
factors  such  as  air  column  stability  or  water  mass  type  can  influence  the 
areal  coverage  of  whitecaps,  photographic  data  was  collectea  over  the  Loop 
Current  from  a  NAVOCEANO  P-3  aircraft.  The  variation  of  whitecap  coverage 
along  a  line  of  closely  spaced  (25  km)  stations  was  determined  and  compared 


DD  ,  :°r7]  1473  EDITION  OF  I  NOV  *5  IS  OBSOLETE 

S  'N  OT  02- L  F -01  4-6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAOe  fWi»n  Dmf  Knltrmd) 


UNCLASSIF 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  ftWiMl  Data  Entarad) 


with  other  aircraft  and  data  buoy  information.  The  data  indicates  that,  on  the 
day  of  the  flight,  white-capping  within  the  boundaries  of  the  Loop  Current 
depended  little  on  the  local  wind,  but  demonstrated  a  noticable  dependence  on 
air  column  stability.  The  strength  and  nature  of  this  dependence  varied  over 
mesoscale  distances.  These  results  indicate  that  microwave  radiometric 
measurements  can  be  sensitive  to  variables  other  than  surface  wind  since 
microwave  brightness  is  quite  sensitive  to  sea  foam.  The  value  of  areal 
whitecap  coverage  as  a  measureable  geophysical  variable  is  also  examined. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PASEfTWi»n  Data  Entarad) 


